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A finite element study is presented of double-diffusive convection driven by g-jitter in a 
microgravity environment. Mathematical formulations are presented and extensive 
simulations are carried out for g-jitter induced fluid flow, temperature distribution, 
and solutal transport in an alloy system under consideration for space flights. Compute 
at ions include the use of idealized single -frequency and multifrequency g-jitter as well 
as the real g-jitter data taken during an actual Space Shuttle flight. Little correlation 
is seen between these velocity components for the g-jitter components studied. The tem- 
perature field is basically undisturbed by convection because of a small Pr number for 
the fluid. The disturbance of the concentration field, however, is pronounced, and the local 
variation of the concentration follows the velocity oscillation in time. It is found that 
although the concentration field varies in both position and time, the local concentration 
gradient remains approximately constant in time. Numerical study further indicates that 
with an increase in g-jitter force ( or amplitude) , the nonlinear convective effects become 
much more obvious , which in turn drastically change the concentration fields. The simu- 
lated results computed using the g-jitter data taken during space flights show that both 
the velocity and concentration become random , following approximat ely the same pattern 
as the g-jitter perturbations. 


INTRODUCTION 

Understanding the effect of gravity perturbation, or g-jitter, on thermally 
induced convection in microgravity is of critical importance in designing and 
optimizing thermal fluids systems for space materials processing. Alloy solidification 
experiments conducted in space vehicles showed that the solute uniformity and 
defects formation in space grown crystals are strongly affected by natural convection 
in the melt pool that arises as a result of the combined action of temperature and 
concentration gradients in the melt and gravity [1—4], These deleterious effects 

Received 29 February 2000; accepted 6 October 2000. 

The authors gratefully acknowledge the support of this work by NASA (Grant No. NCC8-92 and 
NAG8-1693). 

Address correspondence to Professor Ben Q. Li, School of Mechanical and Materials Engineering, 
Washington State University, P.O. Box 642920, Pullman, WA 99164, USA. E-mail: li@mme.wsu.edu 


245 


246 


Y. SHU, B. Q. LI, AND H. C. DE GROH 


NOMENCLATURE 


C, 

heat capacity 

Cffi Cl 

maximum and minimum 
concentration 

AC 

concentration difference 

Co 

dimensionless concentration 

D 

solutal diffusivity 

e 

unit vector 

/ 

dimensionless frequency of the 
acceleration 

r 

frequency of the acceleration 

g 

dimensionless amplitude of the 
acceleration 

s‘ 

amplitude of the acceleration 

go 

earth gravity 

Or s 

solutal Grashof number, 
Gr r = gfijACL'lv 1 

Gr T 

thermal Grashof number, 


Gr s = gpcATU/v 2 

Uj 

unit vectors of the /th, yth 
components 

k 

thermal conductivity 

Li 

length of the square cavity 

P 

discretized nodal pressure array 

Pr 

Prandtl number, Pr-v/a 

Sc 

Schmidt number, Sc = \/D 

T 

dimensionless temperature 

T h . T l 

maximum and minimum 
temperatures 

AT 

temperature between T H and T L 


u 

M, V 

V 

a 

Pc 

Pt 

<t> 

o 

4 

V 

aft 


p 

n 

Subscripts 

ij 

n 

Superscripts 

U 

T 


dimensionless velocity 
x and y component of the 
dimensionless velocity 
volume of the cavity 
dimensionless coordinates 
thermal difTusivity 
solutal expansion coefficient 
thermal expansion coefficient 
shape function for velocity 
shape function for temperature 
and concentration 
shape function for pressure 
dynamic viscosity 
kinematic viscosity 
boundary of computational 
domain 

gradient operator 
density 

computational domain 


/th, yth point 
nth component 


/th, yth component 
matrix transpose 


generally are considered to be a culprit for the nonrepeatability of the quality of 
crystals grown in microgravity [5-8], 

Gravity perturbation or g-jitter originates from a variety of sources including 
crew motions, mechanical vibrations (pumps, motors, excitations of natural fre- 
quencies of spacecraft structures), spacecraft maneuvers and attitude, atmosphere 
drag, and the Earth’s gravity gradient [9], Figure 1 shows the typical signature 
of g-jitter registered by an accelerometer installed on the Space Shuttle during 
an actual space flight. Clearly, g-jitter is random in both time and direction and 
is sharply different from terrestrial environment. Research on g-jitter induced natu- 
ral convection m single component fluid flow systems suggests that convection in 
microgravity is related to the magnitude and frequency of g-jitter and to the align- 
ment of the gravity field with respect to the direction of the temperature gradient 
in melt processing systems [5, 6], Recently, both 2-D and 3-D numerical simulations 
have been conducted to estimate and calculate the adverse effects of time varying 
g-jitter [5, 6, 10-13], The calculations used both idealized single- and multiple- 
frequency g-jitter modulations. These studies showed that the frequency, amplitude 
and spatial orientation of the residual gravity vector all play an important role in 
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determining the convective flow behavior of a melt being processed in space. When 
the residual accelerations oscillate about the positive and negative of an axis, 
the orientation of this direction relative to the density gradient determines whether 
a mean flow is generated in the liquid [5, 6]. Sinusoidal oscillating accelerations 
induce an oscillating convective flow in the liquid. The resultant velocity fields 
oscillate with the same frequency as the affecting gravity field. 

Although these studies are useful in illustrating the basic features of g-jitter 
induced convection in a single component system, very little information seems 
to be available on the fundamental understanding of double-diffusive convection 
in a microgravity environment. An exception is a paper published in a recent con- 
ference proceedings by Bune et al. [14], who discussed double-diffusive convection 
during solidification of HgCdTe. The study of double-diffusive convection can 
be of critical importance in binary alloy solidification systems, because the quality 
of the final products is strongly correlated to the concentration distribution in 
the melt during processing. The main reason for conducting solidification exper- 
iments in space is to eliminate convection induced solutal striation. Several alloy 
solidification experiments are currently being planned for future Space Shuttle flights 
and for the International Space Station under construction. This paper presents a 
numerical analysis of double-diffusive convection induced by g-jitter in one of these 
alloys melt systems in microgravity environment. Toward this end, a simple 2-D 
cavity geometry is considered and the Navier— Stokes equations, energy equations 
and solutal transport equations are solved using the finite element method The 
analyses will be presented of the effects of both the g-jitter data synthesized using 
the Fourier series transformation as well as real g-jitter data taken directly during 
space flights. The information and computation methodology should be of great 
value for gaining an understanding of the basic physics governing the g-jitter induced 
double-diffusive convection, which allows rational guidelines to be developed for 
designing and operating thermal fluids systems for space experiments involving 
binary melt processing. 


PROBLEM STATEMENT 

Figure 2 schematically illustrates the model problem to be studied along with 
the coordinate system for analysis. We consider a square cavity filled with a binary 
liquid metal alloy. Its left side is fixed at a higher concentration C H and right side 
at a higher temperature T H . One of the important characteristics is that the solutal 
Grashof number Gr,y and the thermal Grashof number [Grx be approximately 
the same. The dimension of the square is set to L, as nondimensionalized, equal 
to 1. The upper and lower walls are kept adiabatic. The origin of the coordinate 
system is selected at the lower left corner of the cavity. In microgravity, 
double-diffusive melt convection arises because of the combined effect of the tem- 
perature and concentration gradients and the time varying gravity perturbation. 
The melt flow in turn will alter, in principle, both the thermal and concentration 
fields. The objective of this paper is to develop a basic understanding of g-jitter 
double-diffusive convection and its effect on solutal distribution in the melt mixture, 
through extensive numerical simulations and analyses. 
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Figure 2. Schematic representation of g-jitter driven oscillating convection in a 2-D cavity in microgravity 
environment. 


While g-jitter in general is random in both time and direction, as seen in Figure 
1, in some cases and for a certain period of time, g-jitter may be partially represented 
by a synthesized Fourier series via the fast Fourier transformation, viz., 

gV)=£>;sin(2K/;<*)e (1) 

n= 1 

where g * f J and f* are the amplitude and frequency of the nih acceleration component, 
respectively. Here the subscript * designates a dimensional quality. 

The g-jitter induced melt flow, heat, and mass transport in the system are 
governed by the continuity equation, the Navier-Stokes equations, the energy bal- 
ance equation, and the species conservation equation. For the current system, these 
equations are written in a nondimensionalized form as follows: 


V • u = 0 


( 2 ) 
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du 


— + (u • V)u = -VP + V 2 u - (GttT + GrsC)g(t) 


(3) 


dT 1 , 

— +u*vr =— v 2 t 

at p r 


(4) 


f + u.VC=iv’C 


(5) 


n the above equations, lengths are nondimensionalized by L, velocity by v/L time 
y I, / v , temperature by AT = Th — Tl, and concentration by AC — Cn - Q.. Also 
use has been made of Boussinesq’s approximation, p(T, C) = p Q ( 1 - p T (J - T 0 ) - 
Pc(C ~ C 0 )), and g(f) is nondimensionalized by g 0 = 9.8 m/s 2 . 

To obtain solutions to the above equations, the following boundary conditions 
are applied: 


u — 0 at x — 0 and x = 1 and y — 0 and y = 1 
dT 3C 

~fy = ~fy =0 at ^ = 0and >'=l 
T = 1 C — 0 at x = 0 
r = o C = I at x = 1 


THE FINITE ELEMENT SOLUTION 


TJe above governing equations along with the boundary conditions are solved 
using the standard Galerkin finite element method. The stiffness matrix is obtained 
by using Galerkin’s method of weighted residuals. Following the procedures given 
in [15, 16], the governing equations are recast in the following integral forms 


(1 ^ ‘ V(f>T dV ) U > = " £ (j[ 

(I HT dV ) ^ L+ (l <f>a ‘ V<i>TdV )-(l l V ^ T dV ) P 

+ (I • v<t> T dv ^v, + QV • V0X/ • v<t> T )dv^Uj 

+ (I Grr8i< 1 >dT ') T + (I Gr sg,<i>0 T ^C = j n • T • Ip dS 
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(/„" r dv )% + (/a 9 ”- wr dv ) T+ {ir, ve - veT dv y = -L ,Te ds 

( 8 ) 


(j> r dy )f- + {k“- veT dy ) c+ {Lt e - veTdv ) c = -L qreds 

( 9 ) 

Once the forms of shape functions if>, 0, and ^ are specified, the integrals 
defined in the above equations can be calculated numerically over each individual 
element. Combining the discretized momentum, energy, and solutal transport 
equations gives rise to the following global finite element stiffness matrix equation: 


MOO 
0 N r 0 
0 0 N c 


U 

T 

C 


A(U) + K + l/sEMj'E 7- 


0 

0 


Br 

D r (U) + L r 
0 


Be 

0 



u" 


T 


C 


F 

Gr 
G c 


(14) 


It is noted that in constructing the above element matrix equation, the penalty for- 
mulation has been applied, and P in the momentum equation is substituted by 
(1 /e)Mj'E T U. The assembled global matrix equations are stored in the skyline form 
and solved using the Gaussian elimination method. The above coefficient matrices 
from above are calculated by 


M,= 

f \b\li T dV Nr = 

f ee T dv 


Jn 

Jn 

M= j 

f <p(j> T dV N c = 

f 66 t dV 

J 

n J 

'n 

= f ^r*V0 • V0 r dV 
Jn Sc 

E, = j 

f j ■ V<M r dV L c 

0 

Lr = 

f J_V0.V0 r </F 

In Pr 

A(U) = [ <j)U • V0 r dV 
Jn 

D C (U) 

f- 

> 

S 

II 

D r (U)= f 0u • V0 r dV 

Jn 

Jn 
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= f GT T g<t>6 T dV B c = [ GT S g<t>6 T dV 

JCi Jn 

— ~ I <Jt 6 dS G c = — f qc6 dS 

Jxi 

) = / n • t • i<p dS 
J an 

K (/ = (jf V0 • V< t> T dyS j d “ + jO • V <W • V) dV 


Br 

Gt 

F, 


RESULTS AND DISCUSSION 

The transient finite element computational methodology described above 
enables the prediction of the transient convective flows driven by g-jitter, the tem- 
perature distribution and the solutal redistribution in the melt. The calculations were 
performed using a finite element code developed by our research group [1 6]. The code 
employs the backward Euler (or implicit) time difference scheme, with automatic 
time step control. The computer code was checked against other available commer- 
cial codes for fluid flow and thermal calculations, including FIDAP and FLUENT, 
and analytical solutions for a wide range of flows and thermal problems. Extensive 
numerical simulations were made. A selection of computed results is presented below 
and the parameters used for calculations are given in Table 1. Most transient cal- 
culations are started with an initial thermal and fluid flow field obtained with a con- 
stant gravity field of g = 1 x 10 6 go acting in the direction perpendicular to the 
growth direction to simulate a situation in which a long steady state microgravity 
is perturbed by g-jitter. 3 

Prior to embarking upon numerical simulations, a variety of element types and 
nodal distributions were tested to obtain an appropriate mesh arrangement for com- 
putational accuracy and optimal computing time. The tests included 100, 225 400 
and 625 elements with 9 nodes and 400, 900, 1,600, and 2,500 4-node’ elements! 
For all these cases tested, a spatial distribution of nodes is arranged such that 
the spacing is reduced following the arithmetic progression from the middle of 
the cavity toward the four walls. The conditions used for the testing correspond 
to the maximum possible driving force for double-diffusive convection to be studied 
in this paper. These tests showed that 225 9-node elements with a node spacing 
decreased from the center toward the walls following an arithmetic progressive ratio 
of 0.3 that gave an optimal finite element mesh for the computation in both time and 
accuracy. Any further refinement of the mesh produces the results with an error less 

than 0.2%. The mesh is shown in Figure 3 and is used for all computations presented 
below. 

During space flights, the idealized microgravity level is on the order of 
1 x 10 go and is essentially steady until a disturbance or g-jitter sets in. Thus, 
the transport variable fields may be considered to be in a quasisteady state when 
flights are in perfect condition, thereby enjoying the true benefit of microgravity. 
Figure 4 depicts the quasi steady state temperature, concentration, and fluid flow 
fields in the system under consideration for a typical steady microgravity condition. 
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Table 1. Parameters used in calculations 


Dimensionless parameter 

Value 

Prantal number (Pr) 

0.012 

Grashof number (Grr) 

1.54E6 

Grashof number (Gr$) 

-1.54E6 

Schmidt number (Sc) 

65.6 


Apparently, the convection in the melt is very small, as indicated by the velocity plot, 
and so is the convective transport of energy and species in the system. The velocity 
arises near the hot wall and recirculates downward along the cold wall, forming 
a rotating flow loop. The transport of energy clearly is controlled by a diffusion 
mechanism, which is expected because of the small Prandtl number of the fluid. 
The convection, however, affects the concentration distribution in this case, as is 
evidenced by the distorted isosolutal contour lines. When the system is perturbed 
and g-jitter sets in, the ideal dynamic weightlessness condition is destroyed, and 
the microgravity level is no longer stable (see Figure 1) and changes both spatially 
and temporally, thereby inducing time dependent convective flows and the solutal 
striation. In some cases g-jitter can be represented by a synthesized Fourier series, 
which helps to facilitate simulations and analyses. In what is considered below, 
we present the analyses of the effects on the convection in the double-diffusive system 
that are generated by g-jitter with single-frequency and multifrequency, and by 
measured g-jitter data. 


Figure 3. Finite element mesh for com- 
putation. Nine nodes are used in each element. 
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Figure 4. The steady state fluid flow field and its thermal and solute distribution caused by a constant 
g-jitter (g - 10” 6 go) acting in the v-direction: (a) the velocity profile = 0,012! 1); (b) the isothermal 
profile; and (c) the isoconcentration profile. 


Single-Frequency g-Jitter 

Single-frequency g-jitter represents a somewhat idealized situation in that the 
true g-jitter disturbance consists of a wide range of the frequency spectrum. 
However, a study of this idealized single-frequency disturbance helps us to gain 
a fundamental insight into the basic behavior of the double-diffusive convection. 
Space experiments suggested that melt growth systems are prone to the disturbance 
of the g-jitter components of low frequency. The harmful frequency range lies in 
the range of 0.001 to 1 Hz [5, 6]. Numerical results were obtained for g-jitter induced 
convection and are shown in Figures 5 and 6 for a g-jitter component of 
g = 1 x 10~ 3 gosin(0.27r0* which represents a typical harmful disturbance to a melt 
processing system, with the initial condition given in Figure 4. In Figure 5, the 
snapshots of the oscillating velocity and concentration profiles are given for one 
cycle period while the transient development of field variables is plotted at some 
specific locations in the cavity. As stated above, the temperature profile remains 
approximately the same with very little, unnoticeable oscillation and hence is 
omitted here. Inspection of these results shows that convective flow in the system 
has reached a quasisteady state in a time frame of approximately 3 or 4 periods, 
and afterwards it oscillates with basically the same frequency of the g-jitter 
(dimensionless period = 0.02688). The convection lags in phase, though, the phase 
angle being a function of the amplitude of the g-jitter component, the location, 
and flow field. That the flow oscillates with the same frequency of the driving g-jitter 
may be attributed to the fact that although the main flow is strong, the velocity 
gradient is much smaller, and thus the product of the two contributes negligibly 
to the flow. 

The oscillating flow has a strong effect on concentration distribution. In fact, 
concentration, once it reached a quasisteady state, also oscillate harmonically with 
time with approximately the same frequency over a time-averaged value; again 
the phase angle is a function of the amplitude of the g-jitter component, the location, 
and flow field. For this system, it appears that the concentration oscillation primarily 
sterns from the oscillating velocity field. There is yet no subharmonic (or 
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Figure 5. Evolution of the velocity ( a)-{d) and concentration (e)— (A) fields induced by a single-frequency 
g-jitter driving force in they-direction: (a) t = 14.33163 (C/nvui = 0.6367);, ( b ) t = 14.33163 concentration; 
(i c ) t - 14.33970 (U M - 2.752); (rf) / = 14.33970 concentration; ( e ) / = 14.34507 (U mu = 0.6368); (/) 
/ = 14.34507 concentration; (#) r = 14.35314 (£/„»* = 2.752); and (A) /= 14.35314 concentration. 
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Figure 6. Time variation of velocity components (J, V and concentration with oscillating gravity; (a) (d) 
(g) at (0.07798, 0.07798); (b), ( e ), (h) at (0.9220, 0.5); and (c), (/), (/) at (0.5, 0.9220). 


frequency-doubling) effect. In fact, detailed analyses show that although the con- 
centration field oscillates in time, its gradient remains basically time independent. 
As a result, the product of the velocity and concentration gradient (or the initial 
term) follows the same frequency as the velocity, thereby producing the results 
as shown in Figure 6. This time harmonic oscillation is preceded by a transient 
period, which depends strongly on the location and direction of the concentration 
gradients (see Figure 6). For the points near the vertical wall where the main flow 
is perpendicular to the concentration gradients, the transient period is about 200 
periods. For the points near the horizontal walls where convective flow is parallel 
to concentration gradients, the transient period is as long as 500 periods. The ampli- 
tude of concentration oscillation is also much larger at these locations, which clearly 
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is attributed to the fact that the flow parallel to the concentration gradient is much 
stronger. 

Additional numerical simulations also were carried out by switching off either 
the concentration or temperature field. It is found that the combined 
double-diffusive convection for this particular system is approximately a 
superposition of the concentration and temperature fields, and that both thermal 
and solutal driving forces play almost equally important roles in driving the con- 
vection. 

Extensive numerical experiments also were conducted for various g-jitter con- 
ditions and variations. Results suggest that the evolution and distribution of the 
field variables (velocity and concentration) are sensitive to the initial conditions. 
For example, if the initial fields are generated by a constant gravity field of 
g y = 10~ 3 go> then the change in velocity will take much longer to reach a quasisteady 
state. For the same point near the middle bottom wall, it takes approximately 300 
periods for the velocity to reach the time harmonic quasisteady state. For those 
points near the vertical walls, the transient period lasts about 700 cycles of 
oscillation. 

Detailed analyses of the results indicate that there exists a subharmonic 
oscillation, or the doubling of frequency, in both secondary flow components (or 
velocity component perpendicular to the wall), although very week. The concen- 
tration field, however, does not exhibit subharmonic oscillation until the amplitude 
of g-jitter is increased to 10 _2 go- Under this condition, the concentration contour 
is significantly distorted, and one of such cases is shown in Figure 7. The temperature 
field, however, remains the same until a much larger g-jitter amplitude, which is not 
at all realistic, is applied. The distorted concentration profile, combined with an 
increase in secondary flows, causes a more pronounced nonlinear oscillation (or fre- 
quency doubling) in velocity components perpendicular to the wall; the main flow 
components, however, remain basically the same in an oscillation pattern. Such 
a disturbed flow field causes the concentration distribution to be more strongly dis- 
torted. This is exhibited in both the concentration contour and the concentration 
oscillation at specific locations. As illustrated in Figure 8, the frequency-doubling 
effects, which clearly come from the convective term, start to show up in the velocity 
components perpendicular to the wall and in concentration as well after about 10 
cycles of oscillation. 

Multiple-Frequency g-Jitter 

Let us now turn to the effects of multiple-frequency g-jitter on the convective 
behavior in the system. Toward this end, we select a representative frequency 
spectrum, f — 10“ 3 ~ 1 with \g\ = 10” 4 — lO^go. respectively, which are con- 
sidered to be able to generate detrimental effects in melt processing systems. 
The computed results for velocity and concentration are plotted in Figures 9 
and 10. As expected, the temperature field remains practically unchanged and hence 
is not shown. Detailed examination of the convection patterns shows that the main 
oscillating recirculating flow in the cavity is attributed to the g component with 
the largest amplitude and other components produce only perturbations that are 
superimposed on the main flows. This is further confirmed by the results of velocity 
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( a ) (b) 

Figure 7. The dynamic development of the velocity and concentration fields generated by a 
single-frequency g-jitter driving force (? = 10" 2 g 0 sin(O.2n/)): (a) the velocity profile (f/ m „ = 13,14); 
(b) the isoconcentration profile. 


and concentration variation at specific locations, from which it is clear that the evol- 
ution of the velocity and concentration fields is an approximate Fourier summation 
of velocities generated by each individual g-jitter component. This observation 
should not come as a surprise, in light of the above discussions on results obtained 
for single-frequency components. 


Real o-Jitter Taken from Space Shuttle Flights 

In reality, g-jitter is randomly in nature, as a result of a variety of sources can 
cause the perturbation of the dynamic weightlessness conditions of space vehicles 
in orbit. Numerical simulations were carried out using these real g-jitter data as 
shown in Figure 1 and some of the results are given in Figures 11 and 12. There 
exists no regulated transient pattern or time harmonic motion, and thus the con- 
vective flows in the system must be constantly monitored at any instant. The tem- 
perature field, however, remains the same, again a result of the predominance of 
thermal conduction in fluids. As might have been expected from the single and mul- 
tiple g-jitter analyses, the flow pattern is similar with a single recirculating loop, 
but the time evolution is no longer regular in frequency because of the irregularity 
of the g-jitter driving force typical of the true nature of gravity perturbation in space. 
The time variation of velocities and concentration at specific locations further illus- 
trates that the velocity responds quickly to the change in gravity perturbation, 
whereas the concentration does not. In fact, the concentration still is evolving in 
the transient stage, reflecting the mixed diffusive-convective nature of the concen- 
tration development in the system, even after the gravity perturbation disappears. 
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Figure 8. Time variation of velocity and concentration distribution with single-frequency g-jitter driving 
force = !CT 2 £o sin(0.2*f))' ( a ), (c), (e) velocity components U % V and concentration at (0.9220, 0.5); 
and ( b ), (</), (/) velocity components U, V and concentration at (0.5, 0.9220). 
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(«) (*) 




Figure 9. The dynamic development of the velocity and concentration fields caused by a 
multiple-frequency g-jitter driving force: (a) t = 2.715245 (Umax = 2.709); ( b ) t = 2.715245 concentration 
profile; (c) ( = 2,745353 (On** = 0.5530); and (d) t — 2.745353 concentration profile. 


CONCLUDING REMAKRS 

This paper has presented a finite element study of double diffusive convection 
driven by g-jitter in a microgravity environment. Extensive simulations were carried 
out using both the synthesized single- and multifrequency g-jitter and also the real 
g-jitter data taken during space flights. Results were analyzed for fluid flow, tem- 
perature distribution and solutal transport in an alloy system under consideration 
for Space Shuttle flights. Computed results suggest that for the system studied 
the main flow is in essence a combination of velocity components driven by each 
individual g-jitter component. The flow field is characterized by oscillating 








ionless velocity component V Dimensionless vetoed oo 
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Figure 10. Time variation of velocity and concentration distrubution caused by multiple frequency g-jitter 
driving force: (a t c, e) velocity components U, V and concentration at coordinate (0.9220, 0.5), (fr, d,f) 
velocity components U, V and concentration at coordinate (0.5, 0.9220). 
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Figure 11. The dynamic development of the velocity and concentration fields caused by real g-jitter 
driving force; (a) t = 0.2698170 ( CW = 0. 1 579); (6) / = 0.2698170 isoconcentration profile; (c) 
t — 0.5386340 (t/mu =0.06303); and (d) t = 0.5386340 isoconcentration profile. 


recirculating natural convection, with an oscillating frequency being equal to that of 
the driving g-jitter force. There exists a lag in the phase angle, which is a function of 
position, the local velocity, and the amplitude of the g-jitter force. There exists only 
a very short period of transient time, which lasts about a few periods of oscillation. 
The temperature field is basically undisturbed by convection because of the fact that 
the fluid has a small Pr number. The effect of convection on concentration dis- 
tribution, however, is much more pronounced. The local variation of the concen- 
tration follows approximately the velocity oscillation in time. It is found that 
while the concentration varies in both position and time, the local concentration 
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Figure 12. Time variation of velocity and concentration distribution under real g-jitter driving force: ( a ), 
(c), (e) velocity components U, V and concentration at (0.9220, 0.5); (A), ( d ), (/) velocity components U, V 
and concentration at (0.5, 0.9220). 
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gradient changes remain approximately constant in time. This is combined with the 
local velocity to give rise to the oscillating concentration field with the same fre- 
quency of flow field. The transient period for the concentration field is much longer 
than that for the velocity, however; and can be as long as 500 periods of oscillation 
before the concentration field finally reaches a quasisteady state time harmonic 
oscillation. Numerical results further indicate that with an increase in g-jitter force 
(or amplitude), the nonlinear convective effects become much more obvious, which 
in turn drastically change the concentration fields. The simulated results computed 
using the g-jitter data from real space flights show that both the velocity and con- 
centration are random, following approximately the same pattern as the g-jitter. 
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